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ABSTRACT 

We examine the possibility of observing gravitational lensing in the weak 
deflection regime by the supermassive black hole in the center of the galaxy 
M31. This black hole is significantly more massive than the black hole in the 
center of our Galaxy qualifying itself as a more effective lens. However, it is 
also more distant and the candidate stellar sources appear consequently fainter. 
As potential sources we separately consider stars belonging to the bulge, to the 
disk, to the triple nucleus formed by PH-P2 and by the recently discovered inner 
cluster P3. We calculate the number of simultaneously lensed stars at a given time 
as a function of the threshold magnitude required for the secondary image. For 
observations in the K-band we find 1.4 expected stars having secondary images 
brighter than K = 24 and 182 brighter than K = 30. For observations in the V- 
band we expect 1.3 secondary images brighter than V = 27 and 271 brighter than 
V = 33. The bulge stars have the highest chance to be lensed by the supermassive 
black hole, whereas the disk and the composite nucleus stars contribute by 10% 
each. The typical angular separation of the secondary images from the black hole 
range from 1 mas to 0.1". For each population we also show the distribution of 
the lensed sources as a function of their distance and absolute magnitude, the 
expected angular positions and velocities of the generated secondary images, the 
rate and the typical duration of the lensing events. 

Subject headings: Gravitational lensing — Black hole physics — galaxies: indi- 
vidual (M31) — galaxies: nuclei 



1. Introduction 



As soon as the evidence of a supermassive black hole in the center of the Milky Way 
(identified with the radio source Sgr A*) became overwhelming, the scientific community 
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started to investig ate its ability of acting as a grav itational lens for background stellar sources. 
The first work by I War die fc Yusuf-Zadehl (Il992l ) considered lensing of sources in the imme- 
diate environment of Sgr A* and showed that the background sky should appear slightly 
depleted by a lensing effect. They suggested the possibility that gravitational lensing could 
be strong enough to generate an observable s e conda r y image for several star s at a given 
time. More pr e cise e stimates by IJaroszyhskil (119981 ). [Alexander fc Sternberg] (119991 ). and 
Chaname et al.l (120011 ). point out that about 10 stars belonging to the Galactic bulge should 
give rise to secondary images with < 23 at a given time. Deep and high resolution images 
in the near infrared of the Galactic center are being obtained by several advanced observa- 
tories such as Keck, VLT, and Gemini North. The progress registered by these observations 



encourages to search for signatures o f lensing effects 



two images of a background source (lAlexanderl 1200 



the background stars due to gravitational lensing (INusser &: Broadhurstll2004l ). could help 



3y Sgr A*. For example, identifying the 
), or measuring the astrorn etric shift of 



to determine the position of Sgr A* very accurately. Enhanced microlensing caused by black 



holes surrounding Sgr A* has also been considered (lAlexander fc Loebll2001 



Chaname et al. 



20011 ). A fascinating possibility is offered by stars orbit ing very close to Sgr A*, which have 



now been followed very accurately across s everal years ( 



199' 



Ghez et al. 



Weinberg et al 



199 



20051: 



Eckart et al. 



200 



Eisenhauer et al. 
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Schodel et al. 



20051: 



2002L 



Paumard et al. 
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1997; 



Genzel et al. 



Ghez et 



Reid et al 



al 



2003 



20061 ). In fact, 



2005 



the precise knowledge of the position of the source allows to predict the time, the position 



and the brightness of the secondary image (jPe Paolis et al.l l2003l : iBozza &: Mancinil 12004 
2005h . 



Meanwhile, our kno wledg e abou t the central regions of the Andromeda galaxy (M31) 
has grown considerably. iKentI (119891 ) has provided a detailed brightness and de nsity profile 



for th e bulge and the disk of M31, which is still the reference for present works. iLauer et al. 



( 119931 ) discovered that th e nucleus of M31 is actually constituted by two components, which 
they called PI and P 2. iKing et al.l (1995! ) noticed that P2 is much brighter than PI in 
the ultraviolet. Then Tremain3 (1995) hypothesized that PI and P2 are actually parts of 
the same eccentric disk orbiting a supermassive black hole, with PI being fo rmed by stars 



at the apocenter and P2 being formed by stars at the pericenter. Later on, iBender et al. 



( 120051 ) discovered a bright cluster of young stars (named P3) embedded within P2. They 
also managed to give the most precise estimate for the ma s s of t he central black hole of 
M31, M = lAtoi X 10^ Mq. Recently, IPemarque fc Viranil (12003) have cast doubt on the 
hypothesis that the P3 cluster is made of young stars of spectral class A5 - B5, proposing that 
it is actually made of old stars. On the other hand, IChang et al.l ( 120071 ) have supported the 
young-stars hypothesis suggesting a mechanism for the refuelling of gas in the neighborhood 
of the central black hole. 
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It is very interesting to note that the supermassive black hole in the center of M31 
is nearly two orders of maKnitudes h e avier than the black hole in our Galaxy, estimated 



to 3.6 X 10^ Mq by lEisenhauer et al.l (120051 ) . From the lensing point of view, the greater 
distance of M31 with respect to Sgr A* is therefore partly compensated by the larger mass 
of its central black hole. This motivates a deep investigation of the probability of having 
stellar sources lensed by the supermassive black hole in M31. 

Gravitational lensing effects have been successfully employed in several contexts as a 
tool to study the stru cture of galaxies, large-scale structures and cosmological parameters 



(iKochanek et al.ll2004l ). As we will show in this paper, it is not unlikely that gravitational 
lensing by the supermassive black hole in M31 will follow the wake of its ancestors. With 
a conspicuous number of events it will be possible to undertake a precise reconstruction of 
the mass distribution in the inner core of M31. Gravitational lensing would then provide 
an independent and unbiased method to be crossed with other investigation methods such 
as object counting, spectroscopic or proper motion measurements and so on. Indeed the 
study of the physics of the core of M31 would greatly benefit from these new data, which 
will help to understand the physics of the stellar environment of supermassive black holes 
and possiblly shed light on the true origin of these enigmatic objects. 

In this paper, we shall exploit the present knowledge of the bulge and the disk popula- 
tions in M31 to calculate the expected number of simultaneously lensed sources for a given 
threshold magnitude in either the i^'-band or the \^-band. We will also consider sources 
belonging to the central clusters, modelling PI and P2 as components of a single eccentric 
system and P3 as a separate inner stellar cluster surrounding the central black hole. For 
each of the four populations (disk, bulge, P1+P2, P3), we shall present several probability 
distributions characterizing the properties of the lensing events. The paper is structured as 
follows: in § 2 we indicate our reference models for the source populations and describe their 
features. In § 3 we review some basics on gravitational lensing, with particular reference to 
the black hole in M31. In § 4 we present our estimates for the number of lensed sources at 
any given time for a given threshold magnitude of the secondary image. In § 5 we examine 
the contribution of sources of different magnitude to the total number of events. In § 6 
we give the distribution of events as a function of the source distance. In § 7 we show the 
distribution of the angular positions of the images. In § 8 we present the distribution of 
their apparent angular velocities. In § 9 we estimate the rate and the average duration of 
the lensing events. In § 10 we discuss some issues concerning the identification of the lensing 
events. 5 11 contains the conclusions. 
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2. Source populations 

Our aim is to calculate lensing probabilities for four different populations of stellar 
sources around the supermassive black hole in the center of M31, namely P3, P1+P2, the 
bulge and the disk. We therefore introduce specific models for each population containing 
information on their spatial distribution, kinematic properties and luminosity function, which 
are the necessary elements for the calculation of the lensing probabilities to be treated in 
the succeeding sections. 

For all populations, we choose to normalize the spatial distribution /p^ to unity, so 
that the quantity fp^^x, y, z)dxdydz represents the probability of finding a single star of the 
population Pj in the space element dxdydz. 

We shall present gravitational lensing probability estimates for hypothetic observation 
programs in the l^-band (~ 0.55 /im) or in the i^-band (~ 2.2 /xm). The i^-band presents 
some considerable advantages for this kind of research, since the interstellar extinction is 
lower in the near infrared than in the visible. In fact, in the i^-band the total average 



extinction (M31 + our Galaxy) is Ak = 0.1 (jOlsen et al.l 120061 ). whereas in the V^-band 



we have Ay = 0.31. The latter value can be deduced b y adding t he intrinsic extinction 



of the M31 galaxy Ayint = 0.12 (half the value given by iHaru Il989l . since our sources ar e 
close to the center of M31) to the foreground extinction Ay^ext = 0.19 (ISchlegel et al.lll998l ). 
Moreover, large ground-based telescopes are optimized for interferometry in the infrared 
bands. In the most powerful configuration, they can reach resolutions of order the mas. 
Such capability would be very precious for the detection of secondary images generated in 
gravitational lensing events. Even the forthcoming Jam es Webb Space Telesc ope (JWST) is 



designed to carry instruments for deep infrared imaging JCardner et al.lbooeh . On the other 



hand, the V^-band is less affected by background noise. Furthermore, at the diffraction limit, 
the l^-band enjoys a resolution 4 times better than the i^-band, though large interferometers 
operating in the visible bands are still far to come. 

In any case, by comparing the results in the two bands we get a much deeper under- 
standing on the source selection operated by the gravitational lensing phenomenon. Since 
the stellar populations have different luminosity functions in the two bands, all gravitational 
lensing distributions look different, with some features exalted or depressed. This also allows 
a quick double-check of our results. 

On the basis of this choice, we need the luminosity functions of each source population 
both in the i^-band and in the V^-band. In the i^-band, this function will be expressed 
as np.^Mx), defined so that np.{MK)dMK represents the number of stars belonging to the 
population Pj with absolute magnitude in the range [Mk,Mk + dM^]- The total number 
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of stars belonging to the population Pj is recovered after integration on all magnitudes 



npmK)dMK. 



The same definitions hold for the V^-band, with the obvious changes in the notation. 
The following subsections explain the construction of the luminosity functions and the choice 
of the spatial distribution for each stellar population. 



2.1. P3 



The cluster P3 was discovered by iBender et al.l (120051 ) through spectroscopic observa- 
tions using the Hubble Space Telescope (HST). Although this cluster is embedded within 
P2 (see § 12. 2p . it is characterized by a distinct stellar population and different kinematic 
properties, strongly indicating that it must be considered as a separate entity. 

The characteristics of P3 are consistent with the hypothesis of a circular d isk of stars in 



Keple rian rotation around the central supermassive black hole. In particular, iBender et al. 



J2005h used two mod els to fit the spectroscopic observations: an exponential flat disk and a 
Schwarzschildl ( 119791 ) triaxial model. The best fits were obtained with the flat disk or a thin 
Sc hwarzschild rn o del, w ith axial ratio 0.26. As a reasonable synthesis of the models explored 
by IBender et al.l (120051 ). we choose a classical thick-disk spatial distribution 



/P3(X,|/,Z) 



1 



47rr^3ZP3 



exp 



rv3 



sech 



(2) 



with the disk scale being rp3 = 0.8 pc and thickness = 0.1 pc. The inclination of the 
disk relative to the line of sight is Zp3 = 55°, similar to that of the system P1-I-P2. 

The rotation curve of P3 is symmetric around its center, reaching a rotation velocity 
WP3 = 618 km s~^ and dispersion a-p^ = 674 km s~^. 



For the luminosity function of P3 we have followed IBender et al.l (l2005r). who use the syn- 



thetic color-magnitude diagram generated by the program lAC-STAR by lAparicio fc Gallart 



( I2OO4J ). We have run this program with the same parameters, accepting the hypothesis that 
P3 was generated by a single starburst that o ccurred 200 Myr ago in a gas cloud with 
solar metallic i ty (se e iDemarque fc Viranil |2007| for an alternative proposal). We adopt a 
Kroupa et al.l (Il993l ) initial mass function. The algorithm has returned us 10^ stars with 
Mk < 6. Such a number turns out to be insufficient for a substantial covering of the giant 
stars branches. We have thus generated a second sample of 10^ stars with Mk < 1 and 
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combined the two samples with appropriate weights. The two samples are shown together 
in two different color-magnitude diagrams in Figure [H The cut-off at = 6 does not 
translate into a sharp cut-off in the l^-band. However, the sample is completely unaffected 
for My < 8, which we take as our cut-off in the luminosity function in the ^-band. It is 
important to stress that gravitational-lensing observations select very bright sources with 
much higher probability than faint ones. For this reason, we do not pay much attention 
to the completeness of our simulated samples of stars on the low-luminosity side. In the 
same way, issues concerning the lower cutoff of the initial mass function are not relevant for 
us. All gravitational lensing probabilities are practically insensitive to changes in the source 
distribution function at low luminosities. This statement is also supported a posteriori by 
the distributions presented in § |5l 

Note that the red giant branch (RGB) lies at almost fixed magnitude in the l^-band, 
whereas it spans 6 magnitudes in the K-haxid up to the tip. 

The normalization of our sample of P3 stars has been performed introducing a factor 
N multiplying the total number of stars. The value of has been deduced by comparing 
the total m agnitude in the \^-b and of our artificial set with the total magnitude of P3 as 



deduced by iBender et al.l (120051 ) . being My^ps = —5.7 



Mv,P3 = -2.5 logio N - 2-5 logio Afv,P3-sample, (3) 

where My^P3_sampie is the unnormalized absolute magnitude of our combined sample. Also 
this normalization procedure on the number of stars deduced by the total luminosity is 
largely insensitive to the abundance of low- luminosity stars. 

Finally, the same set of stars has been used to build a binned luminosity function in the 
A'-band and in the V^-band, normalized in the way just described. The values of these binned 
functions are shown in Tables [1] and [21 together with those of the other source populations, 
to be discussed in the succeeding sections. 

The different orientation of the RGB in the V^-band and in the A'-band is reflected in the 
luminosity function, which is zero in the l^-band up to My — —2.75, whereas it is already 
non-zero in the A'-band at Mk — —9.75. The A'-band luminosity function clearly shows a 
red clump at Mx — —3.75, which is not as clear in the l^-band, where it overlaps the stars at 
the turn-off point (TOP) of the main sequence (MS). It is worth noting that the real number 
of stars in P3 should be of order a few thousands, whereas we have simulated 2 x 10^ stars 
in order to have a statistically significant sample. 
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Fig. 1. — (a) (Mr-, My — Mk) color-magnitude diagram for P3; the dashed hne represents 
the boundary between the two samples of stars re-combined with different weights in the 
luminosity function, (b) (My, My — Mj) color- magnitude diagram for P3. 
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Table 1: Binned luminosity functions in the ii'-band for P3, P1-I-P2, the bulge and the disk of 
M31. As each bin spans 0.5 magnitudes, the number of stars in each bin is just the tabulated 
value of Up. multiplied by 0.5. 
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Table 2: Binned luminosity functions in the F-band for P3, P1-I-P2, the bulge and the disk 
of M31. 
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2.2. P1+P2 



On a scale a bit larger than P3, the supermassive black hole in M31 is surrounded by 
two larger clusters, called PI and P2. Actually, PI and P2 are parts of the same eccentric 
disk orbiting the central bla ck hole, with P2 being composed of star s at the periapsis and 
PI by stars at the apoapsis (jTremaindll995l : iPeiris &: Tremainell2003l ). 



The spatial distribution of the whole system was obtained by iPeiris fc Tremaind (120031 ) 
fitting the brightness profile of a simulated sample of 10'' stars with orbital parameters 
randomly chosen from some suitable distributions. We have followed their steps, simulating 
the same number of stars using the best fit parameters of their non-aligned model, which 
provides an excellent fit to the observed brightness profile of P1+P2. In Figure [2^ we show 
the brightness profile t hus obtained. For a l l the details relative to the eccentric disk model, 
the reader is referred to lPeiris fc Tremaind (120031 ). To give an idea of the shape of the cluster 
P1+P2, we just mention that the scale of its extension is fixed by the parameter = 1.37 
pc, its inclination along the line of sight is 6i = 54.1°, its thickness is roughly controlled by 
the combination aocr^, with = 24.6°. For every star, the eccentricity is a function of the 
semiaxis of the orbit with several parameters. It is roughly peaked at e = 0.5. 

Binning our simulated distribution of stars, we then easily obtain the spatial distribution 
function fpi^p2{x,y,z), to be used in the lensing calculations. 

In the kinematic calculations, since the stars behind the black hole are closer to P2, we 
have used the aver age rotation velocity of P 2, estimated as Vp2 = 220 km s~^ with dispersion 
(Tp2 = 100 km s~^ ( Peiris fc Tremaind l2003l ). 



The characteristics of the stars composing P2 and PI are quite similar to those of 
typical bulge stars. We have therefore used the bulge sample illustrated in § 12. 3[ Here we 
j ust mention that we hav e normalized this sample in the V^-band, using Vpi+P2 = 12.55 mag 
( jPeiris fc Tremaindl2003l ). With this normalization, we have built the luminosity functions in 
the i^-band and in the l^-band, as shown in Tabled! Some considerations on these functions 
are included in the following section. 



2.3. Bulge 



The bulge of M31 has been studied in great detail by iKentI (119891 ) . who traced precise 
luminosity and density contours. By interpolating these contours, it is possible to build a very 
accurate spatial distribution. We shall refer to this distribution as /suige- The inclination of 



the plane of symmetry of the bulge is ib 



ulge 



77°. The bulge is assumed to have negligible 



-4 



-2 



2 
X (pc) 



4 




Fig. 2. — (a) Reconstructe d briKhtness profil e of th e central cluster PH-P2, obtained fol- 
lowing the prescriptions of iPeiris fc Tremaind (120031 ): north is up and east is left; the cross 
indicates the position of the central black hole. PI is the brighter cluster, whereas P2 is the 
smaller one appearing closer to the black hole, (b) Spatial distribution of the cluster PH-P2 
along the line connecting the observer with the black hole. The positive Z-axis points toward 
the observer. 
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rotation and constant dispersion velocity (XBuige = 160 km s ^. 

In order to gen erate viable luminos i ty fun ctions for the bulge, we have used the pro- 
gram lAC-STAR by lAparicio fc GallartI (120041 ) in this case as well. In particular, following 
Sarajedini fc Jablonka ( 20051 ) . we have supposed that the bulge (and also PH-P2) were gen- 
erated by a single starburst occurred 12.6 Gyr ago with metallicity following a closed box 
law 

dZ / ' 

with Zq = and yield y equal to the solar metallicity Zq = 0.019. With these specifications, 
we have generated 10^ stars with Mk < 6 and 10^ stars with Mk < 1. Before combining the 
two samples, we have randomly removed stars a t the lower end of the ni e tallici ty distribution 
{Z < 0.004), in order to reproduce the data of ISarajedini fc Jablonkal (120051 ) as accurately 
as possible. This additional cut has reduced the two samples by roughly 10%. Figure [3] 
shows two color-magnitude di agrams of the combine d samp le of bulge stars. By comparing 
with the real ones studied by ISarajedini &: Jablonkal ( 20051 ). it is possible to appreciate the 



accuracy of our sample. 



The luminosity profile introduced by lKent (119891) is obtained in the r-band of the Thuan 
and Gunn filter set. Assuming r — K = 2.9 ( lOlsen et al.ll2006l ). we can directly normalize 
the luminosity function in the i^'-band and consequently in the l^-band. The final results 
are included in Tables [1] and [2j 

The bulge RGB lies in the range -I > Mk > -7.5. The TOP is at Mk ^ -3 rather 
than Mk = as in P3. The red clump is peaked at Mk = —1.25 and is also evident in the 
V^-band at My = 0.75. 



2.4. Disk 



The disk spatial distribution is modelled by 



fDisk{x,y,z) 



1 



47rrf)i3k^Disk 



exp 



f Disk 



sech 



^Disk 



(5) 



with disk scale rnisk = 5 .5 kpc, thickness zoisk = 0.3 kpc and inclination ZDisk = 77° 
dWidrow fc Dubinski|[2005l ). 

The r otation velo city of the stars in the disk is fDisk = 250 km s~^ with a negligible 
dispersion (lKentlll989l ). 
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L i i 1 i 



Fig. 3. — (a) {Mk, Mv — Mk) color-magnitude diagram for the bulge and the cluster P1-I-P2; 
the dashed line represents the boundary between the two samples of stars re-combined with 
different weights in the luminosity function, (b) (My, My — Mj) color- magnitude diagram 
for the bulge and P1-I-P2. 
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The disk star formation history is quite different from the one f ollowed by the b ulge 
and presumably by th e central cluster stars. Following iHodgd ( 119891 ) . IWilliamj (120021 ) and 
Bellazzini et al. J2003h . we take a simple Population I model with a constant star formation 
rate throughout the history of M31 and constant solar metallicity. With these assumptions, 
we have run the program lAC-STAR obtaining 10^ stars with Mk < Q and more 10^ stars 
with Mk < 1. The color-magnitude diagrams so-obtained are shown in Figure |H 

The normalization of our sample of stars has been deter mined comparing th e tota l 
luminosity in the i?-band to the total luminosity as indicated by lWidrow fc Dubinskil (120051 ) 
-Roisk = —21.4. Finally, the correctly normalized luminosity function in the i^-band and in 
the l^-band have been included in Tables [1] and [2l 

With respect to the other populations, the disk has a complete main sequence extending 
up to Mk — —7 and My ~ —7. Moreover, there are also some Red Supergiants (RSG) 
at Mk = —10 and My = —8. The existence of these supergiants mainly affects the V- 
band luminosity function, since it makes the disk luminosity function start at much brighter 
magnitudes with respect to the other populations. The red clump is peaked at Mk = —1.75 
and My = 0.75, though almost contiguous to the MS in the V^-band. 



3. Basics of gravitational lensing 

In this section we shall briefly recall some definitions used in standard gravitational 
lensing. This will prepare the ground for the analysis to be presented in the following 
sections. 

The supermassive black hole in the center of M31 can be modelled as a simple point 
l ens, with mass M = 1.4 x 10^ M©, placed at a distance Dql = 760 kpc from the Sun 



(iBender et al.ll2005l ). We define the optical axis as the line joining the observer to the lens. 



For a source behind the black hole at distance Dls from it, the lens equation takes the form 

/? = ^-f, (6) 

where f3 is the angle between the line joining the observer to the source and the optical axis, 
6 is the angle formed by the observed image with the optical axis and 



V Doi^Dos 

is the Einstein angle and Dqs = Dql + -Dls is the distance from the observer to the source. 
As the source distances considered in this paper are at most of the order of a few kpc behind 
the lens, we generally have Dqs — Dql- 
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Fig. 4. — (a) {Mk, My — Mk) color-magnitude diagram for the disk; the dashed hne repre- 
sents the boundary between the two samples of stars re-combined with different weights in 
the luminosity function, (b) {My, My — Mj) color-magnitude diagram for the disk. 
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Fig. 5. — Magnification of tlie secondary image as a function of the source position normalized 
to tlie Einstein angle. 



Solving the lens equation for 6*, we obtain the position of the two images 

1 



These two images are magnified by a factor 

m2 + 2 



^ 1 



(9) 



where u = (3/9-E,is the source angular position normalized to the Einstein angle. 



Figure shows the magnification of the secondary image as a function of the normalized 
source position. It clearly shows how the dependence of the magnification on u changes from 

to as we go from the regime of good alignment {u <^ 1) to the regime of bad 
alignment (m ^ 1), the Einstein angle representing the scale of the transition between the 
two regimes. 

When the source is far from the optical axis {(5 ^ ^e) the secondary image has very 
low magnification and becomes unobservable. Our aim is to determine the probability of 
seeing the secondary images of sources in M31, which will typically become observable when 
j3 is of the same order of magnitude as 9e- We then note that the angular separation of 
the secondary image from the central black hole is always of the order of the Einstein angle 
for /? ^ 6'e. The order of magnitude of the Einstein angle in the physical situation we are 
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considering is 

The Einstein angle is of order 0.1" for sources as far as a few kpc, whereas it drops to a 
few mas when the source is closer than 100 pc. We can thus conclude that in order to observe 
the secondary images of sources closer than 1 kpc to the central black hole we need very high 
angular resolution, only achievable by interferometric techniques, such as those employed in 
the VLTI, Keck and LBT. This gives an advantage to source populations distributed on a 
larger scale (bulge and disk) with respect to the central clusters P1+P2 and P3. However, 
distinguishing a secondary image at 0.1" from the central black hole from other sources is 
by no means easy and certainly demands high resolution as well. 

Since we are speaking about lensing by a supermassive black hole, it is interesting to 
check that the weak deflection paradigm holds for all the interesting events. Higher order 



images formed by photons performing loops around the central black hole show up 



close to the black hole shadow bord er, which has a radius of 6sh = 9 /ias (jParwinl Il959 



very 



Virbhadra &: EUid l2000l : iBozzal l2002l ) . In order to detect such images, a much greater effort 
with respect to ordinary lensing images would be necessary. For this reason, we will not 
consider them in the present analysis. 

It remains to check whether the secondary image can always be correctly described in 
the weak deflection limit. This can be verified by comparing the position of the image 6- 
with the radius of the shadow border 9sh- Approximating 6^ by 9e, the ratio OsIi/Oe gives an 
estima te of the error we commi t by neglecting the next to leading order term in the deflection 



angle (IKeeton fc Fetters! l2005l ). It is easy to calculate that the error is 2% at Dls = 0.1 pc 
and 6% at Dls = 0.01 pc. As it will be evident in §[21 where we calculate the distribution of 
the events as a function of the source distance, only the low Dls tail of the F3 distribution is 
affected by these errors and then only marginally. Considering the low relevance of this tail 
in the total estimates and the uncertainties in the F3-population modelling, we will simply 
ignore any strong deflection effect in our analysis. 



4. Number of lensed sources 

This section contains the main calculation of this work, namely the number of expected 
lensing events at a given time for each source population around the supermassive black hole 
in M31. We first introduce the methodology that we have followed and then present our 
estimates for each population at the end of the section. We will focus on the i^-band, with 
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analogous considerations holding for the \^-band. 

Consider a source of absolute magnitude Mk- The observed magnitude K is 

ir = M;, + 51ogioT^ + A^, (11) 
lU pc 



where = 0.1 is the extinction in the iiT-band (jOlsen et al.ll2006l ). If the source suffers 



gravitational lensing, its secondary image has apparent magnitude 

K_ = ir-2.51ogio/i-, (12) 

with given by equation ([9]). 

If we fix a threshold magnitude Kthr for the detection of the secondary image, only the 
sources sufficiently magnified will have a secondary image with K_ < -ft'thr- The minimum 
magnification needed to bring a source with absolute magnitude Mk above threshold can be 
found by simply inverting equation f|T2l) with K_ = Kthr 

where AK = Kthr — Mk- 

As the magnification is a function of the normalized source position angle u through 
equation ([9]), there exists a limiting value for u (which we shall indicate by uz) such that a 
source with absolute magnitude Mk has a secondary image just at the threshold value Kthv 
This can be found inverting equation IQ 

uz(/i-) = [2/i_(l + /i_) + (1 + 2/i_)v//i_(l+/i_)] . (14) 

Defining j3z = O^uz, all sources with P < Pz have a secondary image brighter than -ft'thr- 
Finally, we can also define the radius of the lensing zone at distance -Dls, as the radius of 
the circle containing the sources with magnitude Mk whose secondary image is magnified 
above -ft'thr- Of course, the radius of this circle is simply Rz = PzDqs- For any value of 
the distance -Dls? only sources with r < Rz have an observable secondary image. As Z^ls 
varies, we can thus define a lensing zone with radius i?z(-DLs) centered on the optical axis 
and that contains all the sources with absolute magnitude Mk that give rise to observable 
gravitational lensing effects. More explicitly, the radius of the lensing zone is 

i?z = 0.16pc -i^ Mz(/i-(Air)), (15) 



kpc 

which clearly shows that the radius of the lensing zone grows with the square root of the 
distance of the source from the black hole. The dependence on the absolute magnitude of 
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the source and the threshold fixed by observations is stored in the function ui {fi^{AK)). 
The brighter the source and the fainter the threshold, the larger the lensing zone. 

The probability of finding a single star of absolute magnitude in the lensing zone 
defined by the threshold magnitude i^thr is obtained integrating the spatial distribution of 
the population to which the source belongs in the domain contained within the lensing zone 

oo Rz 2tt 

np.{MK,Kti,r) = j dDi^s jdrr j rf^/p^ pLS, r, 0), (16) 



where we recall that Rz is a function of -Dls? Mk, and -ft'thr- In the evaluation of the spatial 
distribution one must take care of the correct geometric orientation in the space of the 
population considered, as specified in § [2l 

Finally, the total number of lensed sources with a secondary above threshold is 

iVp, (i^thr) = j np^ iMK)np^ {Mk, Kttr) dMK. (17) 

Figure shows the estimated number of lensing events for each of the four source 
populations considered in this paper as a function of the threshold magnitude i^thr- Repeat- 
ing the same steps, we can get a similar plot for observations lead in the l^-band, with a 
threshold magnitude Vthr as shown in Figure [6)d. 

We see that P3 gives a practically negligible contribution. The stars in P3 are too few 
and too close to the supermassive black hole. The lensing zone is indeed too restricted at 
distances of the order of a pc to get a sizeable number of events. 

The situation for P1+P2 is much better, because this cluster extends to larger radii 
and is more populated. However, we have to push the threshold magnitude to i^thr = 27 or 
Vthr = 30 in order to have at least one expected event. The estimated number of events for 
the disk is very similar to that of P1+P2. 

The bulge is by far the best reservoir of good sources for gravitational lensing by the 
central black hole. The number of expected events is already larger than one at K^ia = 24 
or Vthr = 27, reaching more than one hundred at K^hv = 30 or Vthr = 32. Comparing the 
plots in the i^-band and the V-band, we can note that redder populations such as P1+P2 
and the bulge are slightly depressed when going from the i^-band to the V^-band. The disk 
supergiants, not present in other populations, keep the the disk number of events higher at 
low values of V^hr with respect to the other populations. 

Roughly, the plot in the V-band is very similar to the plot in the i^-band, but is shifted 
to higher values of Vthr by 2.5 magnitudes. This apparent gap between the two bands is 
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Fig. 6. — (a) Number of expected lensing events as a function of the threshold magnitude 
for the secondary image in the K-hand. (b) Number of expected lensing events as a function 
of the threshold magnitude for the secondary image in the V^-band. Note the different range 
in the ordinate axis in the two plots. 
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actually completely recovered when one compares the signal to noise ratio (SNR) in the 
two bands. In fact, a s ource with V — K = emits 5.77 photons in the \^-band for each 



photon in the i^-band (jCojg |200l|) . In the background-limited regime, the expected noise 
due to the 4 populations considered in this paper can be easily estimated thanks to the 
spatial distributions and luminosity functions introduced in § [2l Indeed we find that for 
images with the same FWHM the SNR in the V"-band is 2.5 magnitudes better than in the 
i^-band. This justifies our choice to plot our expectations in the intervals 20 < i^thr < 30 
and 23 < Vthr < 33. Furthermore, if the diffraction limit in the two bands is reached with the 
same aperture, one gets an additional bonus of 1.5 magnitudes for the ^-band. However, at 
the present time, the largest interferometers such as VLTI, Keck and LET are not designed 
for observations in the l^-band and therefore only the i^'-band can take advantage of the 
resolutions available at such long-baseline facilities. 

The expectations plotted in Figure [6] are also summarized in Tables [3] and HI Recall 
that these estimates give the number of secondary images above threshold simultaneously 
present at a given time. See §[8] and §[9] for a discussion on their evolution in time. It should 
be noted that the estimates presented in this section do not take into account the angular 
separation of the secondary images from the central black hole. § [7| is devoted to that issue. 



5. Magnitude of the lensed sources 

The number of expected events as calculated in the previous section is obtained by 
integrating over the whole range of possible magnitudes for the sources. It is interesting 
to evaluate the contributions of sources with different magnitudes to the final result. In 
Figure [7| we show the integrand of equation (fT7|l for all four populations as a function of 
the source magnitude Mk and the threshold magnitude Kthr- For each value of Kthr the 
function dNp. / (IMk has been normalized to unity. Figure [8] shows the same distributions in 
the l^-band. 

Schematically, we can say that low luminosity stars are more numerous, whereas brighter 
stars enjoy a larger lensing zone and consequently a larger lensing probability. The surfaces 
in Figures [7] and [8] are the outcome of the interplay of these two opposite tensions. At lower 
thresholds, the distributions are peaked on the most luminous stars. As we increase the 
threshold, less luminous populations become predominant. 

As regards the f^-band, in the distributions shown in Figure [7] we can clearly identify 
the tip of the RGB in the high luminosity peak dominating at lower values of the threshold. 
At intermediate values of -ftTthr, the red clump gives the highest contribution. At fainter 
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Table 3: Expected number of events for sources in the bulge, the disk and P1+P2 for different 
threshold magnitudes in the K-hand for the secondary image. 
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Table 4: Expected number of events for sources in the bulge, the disk and P1+P2 for different 
threshold magnitudes in the y-band for the secondary image. 




Fig. 7. — Relative contribution of sources with different absolute magnitudes to the total 
number of expected events for threshold magnitudes ranging from 20 to 30 in the K-hemd. 
(a) P3; (b) P1+P2; (c) bulge; (d) disk. 
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Fig. 8. — The same as Figure [7] in the IZ-band. 
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thresholds, the stars at the TOP become significant and finally overtake the red clump. 
It is interesting to note that the TOP peak gradually shifts to dimmer values of Mk as 
K^Ya is increased and less luminous stars of the main sequence come into play. In the disk 
distribution, the brighter main sequence stars and the RSG also contribute to the first peak 
at low threshold values. 

In the V^-band, the difference between the four populations is made more manifest. We 
can note that P1+P2 and the bulge are characterized by just two peaks: the one at higher 
luminosities is due to the red clump at My — and the dimmer at My ~ 4 is made up 
of stars at the TOP. The RGB is almost horizontal and therefore there is no other peak at 
lower values of My. In P3 the red clump and the TOP give rise to two consecutive peaks 
at My ^ —1.5 and My ~ respectively. As the threshold is increased, the distribution 
fiattens because of the contribution of progressively fainter stars coming into play. The disk 
distributions start with the peak at My ~ —8 due to the RSG, which are absent in the other 
populations. These stars give rise to the different behavior of the disk at low thresholds, 
already discussed in the previous section. At intermediate thresholds the red clump becomes 
dominant and is later followed by the TOP contribution, which first forms a shoulder to the 
red clump peak and then becomes dominant. 

Another important point can be anticipated noting that the distributions for P1+P2 
and the bulge are always practically identical. This proves that differences in the spatial 
distribution have very little influence in selecting the class of sources (luminous or dim) for 
gravitational lensing. This point will be better explained in the next section, where we show 
that the magnitude and distance distributions can be factorized in a flrst approximation. 

As a flnal consideration, we see that the contribution of stars at Mk ^ 6 or My ~ 8 
is always negligible for the thresholds considered in this paper. This provides an a poste- 
riori justiflcation of the irrelevance of the faint end issue of the luminosity function in our 
calculations, as anticipated in § [2j 

6. Distance of the lensed sources 

The four source populations considered in this paper flU different regions of M31, at 
various distances from the center. It is interesting to calculate the distribution of the sources 
suffering gravitational lensing as a function of their distance from the central black hole -Dls- 
Actually, since the distributions are spread over several orders of magnitudes in distance, 
we prefer to present the distributions in log^^Q Dls- In order to get this distribution for each 
population, we return to equations ffTBl) and ffTTj) and reverse the order of integration in Z^ls 



- 26 - 




Fig. 9. — Distributions for the distance of the lensed sources from the black hole. 



and Mk 

dNp. 



oc Z^LS j dMK np^iMj,) j dr r j d<Pfp^{D^s.r,^). (18) 



The distributions so-obtained are shown in Figure [HI after having been normalized to 
unity. The plots are obtained for fCthr = 24, but the distributions for different thresholds 
and even for different bands are practically indistinguishable. This is a consequence of the 
fact that the radius of the lensing zone Rz is always much smaller than the scale of variation 
of all spatial distributions. Then it is possible to approximate the spatial distributions in 
the integrals fll6p and fllSp with the values they assume at the center of the disk of radius 
Rz and replace the integral over r and by TrR^fp. (Z^ls? 0, 0). Since the dependence on Dls 
in the radius of the lensing zone factors out, being always V-Dls? the spatial distribution 
and the magnitude distributions are effectively factorized. This decoupling is not possible 
when we consider the angular position and the velocity of the images, as will be shown in 
the succeeding sections. 

The maximum lensing probability is reached by P3 at Z^ls = 0.19 pc, by P1+P2 at 
Z^LS = 3.7 pc, by the bulge at Z^ls = 0.48 kpc and by the disk at Z^ls = 1-5 kpc. With 
respect to the other populations, the distribution of the bulge events is less localized and 
presents a first subpeak at 5 pc. This is a consequence of the fact that although the bulge is 
sharply peaked at the center of the galaxy, it has a long tail decreasing at a lower rate with 
respect to the exponential tail characterizing the other distributions. Moreover, the other 
populations are more or less flattened on planes that never contain the line of sight, whereas 
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the bulge has a more spheroidal shape. 



7. Position of the images 

From the observational point of view, it is very important to determine the expected 
angular distance from the central black hole of the secondary images we are looking for. 
In practice, the choice of the facilities to employ in the search for these images is heavily 
influenced by the resolution required to resolve them. In this section, we derive the distribu- 
tion of the secondary images generated by gravitational lensing with respect to their angular 
distance from the central black hole. 

Suppose we want to calculate the fraction of sources belonging to a given population 
that generate a secondary image at an angular distance in the range [9, 9 + d9]. The angular 
position of these sources is determined through equation ([H]) as a function of 9 and Dls- 
We thus have to sum up the contributions of sources at different -Dls- Considering that 
\13{9, -Dls) I is a growing function of -Dls, the first contribution comes from a source at distance 
Drmn{9) such that 1/9(6', -Drnin)! = 0. This distance can be exphcitly calculated using equations 
([6]) and ([7j). The last contribution comes from sources at the border of the lensing zone, 
whose distance D^g_^{9) is such that |/3(6', -Dmax) I = /^zl^i^:, -^thr, -Dmax)- This distance is 
also easily calculable (note that the absolute value is necessary because we are considering 
secondary images, i.e. images with 9 < 9e)- Finally, for each source distance in the range 
[-Dmin, -Dmax], wc must sum up the contributions of all space elements lying on a circle of 
radius -Dqs 1/5(6*, -Dls) I centered on the optical axis at distance -Dls from the black hole. 

As in the previous section, we prefer to calculate the distribution in log^Q 9. Taking 
proper account of all the Jacobians, we have 

ax ^ 

cx9 j dM^np^iMK) j dD^s DlM(^,D^s)\ (^2 + ^ j d<Pfp^{D^s, P, 4 



dlogio^ 



(19) 



Equation (fT9|) is for observations in the K-hand. The distribution for observations in 
the V^-band can be calculated similarly with the obvious replacements. Note that now it is 
impossible to factor the dependence on -Dls, -^thr and as done in the previous section. 
As a consequence, the distribution of the angular positions of the images strongly depends 
on the threshold magnitude chosen. 

In Figures [10] and [11] we show the distributions (normalized to unity at each threshold) 
for each source population and for values of the threshold magnitude ranging from 20 to 30 in 
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Fig. 10. — Distributions of the angular distances of the gravitational lensing secondary 
images from the black hole for the four source populations at different values of the threshold 
magnitude in the i^-band. (a) P3; (b) P1+P2; (c) bulge; (d) disk. 
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11. — The same as Figure fTOl in the \^-band. 
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the i^'-band and from 23 to 33 in the \^-band respectively. The images of sources belonging 
to P3 have angular distance in the range [0.0001", 0.001"], those belonging to P1+P2 appear 
at a few mas from the central black hole, those in the disk have angular distance in the 
range [0.01", 0.1"], whereas sources in the bulge can generate images with a larger spread of 
angular distances, from 0.001" to 0.1". In general, as we increase the threshold magnitude, 
the distributions move to shghtly lower values of 9. This can be understood by the fact 
that lowering the threshold we accept fainter images generated by sources more distant from 
the optical axis, whose secondary images appear closer to the black hole. In practice, as 
a general rule, we learn that, in order to catch all fainter images one gets by performing 
deeper observations, one also needs better resolution. This shift is partly compensated by 
the existence of numerous intrinsically dim stars, which anyway need a very good alignment 
in order to generate a visible secondary image. These dim stars tend to keep the distribution 
at the highest possible value of 6', essentially fixed by the Einstein angle at the typical 
distance scale of the population. 

Apart from this general behavior, we can clearly identify some of the populations dis- 
cussed in §[5l For example, in the distribution for P3 in the /T-band we can clearly distinguish 
three peaks at high threshold, corresponding to the tip of the RGB, the red clump and the 
TOP from left to right. These peaks are more smoothed in the distributions for P1+P2 and 
the disk. In the bulge, the double peak structure is mainly determined by its spatial distri- 
bution, which presents a long tail extending up to a few kpc. This long tail is responsible for 
the peak at ^ ~ 0.01". This can be seen by the fact that at i^'thr = 20, at which we expect 
a single peak due to the tip of the RGB, we already have two evident peaks. Moreover, 
the same structure can be seen in the V^-band distribution, notwithstanding the different 
magnitude distribution. 

The V^-band distributions appear to be less structured with just a very shght shift to 
lower values of 9 with the increase of Vthr- 

Summing up, we note that imaging the central regions of the nucleus of M31 with a 
resolution of the order of a few mas, as could be possible with the Keck or the LBT, will 
allow to catch secondary images of sources in the disk and a large part of those in the bulge, 
provided one reaches a faint enough Kt\a (see discussion in § [TU]) . 

8. Velocity of the images 

Up to now, we have just calculated the number of gravitational lensing events that 
are observable at a given time. It is very important to understand the timescale of the 
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variations of the geometric configuration of such lensing events. In fact, the strategies for 
the observations tightly depend on these timescales and can be very different for static or 
nearly static configurations as opposed to events with any secular development. As will be 
discussed in § [101 the detection of the motion of the images might be of key importance to 
identify genuine gravitational lensing events. This section is devoted to the presentation of 
the distributions of the expected velocities of the secondary images below a given threshold 
magnitude. 

These distributions are obtained by a Montecarlo procedure. For each value of Mk 
and Kthr, we have generated 10^ stars, randomly choosing their position within the lensing 
zone corresponding to these magnitudes. The velocities of these sources have been randomly 
extracted from a two-dimensional gaussian distribution centered on the rotation velocity vp. 
with dispersion ap.. The values of these parameters are specified in §[2] for each population. 

Once the sources have been generated, we have calculated the apparent proper mo- 
tions of the corresponding secondary images by differentiating the standard formula for the 
position of the secondary image. The result is 



where ve is the modulus of the time derivative of the angular position of the secondary 
image, v is the modulus of the transverse velocity of the source and is the angle between 
the velocity and position vectors of the source. 

The image velocity decreases as when the source is far away from the optical axis. In 
this limit, the secondary image is also very faint and close to the black hole. In the opposite 
limit, when m — > 0, the velocity diverges (except for the case = 0). This is very well-known 
in microlensing studies, as the secondary image moves very rapidly when the source is at the 
closest approach distance. 

With the velocities of the images thus calculated, we have constructed the distributions 
shown in Figures [12] and [13] for the i^-band and V^-band respectively. The proper motion 
has been expressed in milliarcseconds per year. 

The expected apparent proper motions of the secondary images cover a very large range. 
All four populations generate images whose velocities generally lie between 0.01 and 0.1 
mas/y, with very long tails extending one order of magnitude above and below these limits. 

Going from lower to higher thresholds, all distributions significantly drift to lower values 
of the proper motion. This is clear from the fact that, at bright thresholds, very good 




(20) 
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Fig. 12. — Distributions of the velocities of the gravitational lensing secondary images from 
the black hole for the four source populations at different values of the threshold magnitude 
in the K-hand. (a) P3; (b) P1+P2; (c) bulge; (d) disk. 
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Fig. 13. — The same as Figure [12] in the y-band. 
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alignments generate higher velocities as it is evident from equation fl20|) and viceversa. Then, 
going to fainter thresholds, the lensing zone is enlarged, and less aligned sources participate 
in the distribution pulling it to lower values of vq. This is particularly evident in the \^-band. 

Thanks to the stronger dependence of the velocity distributions on the threshold mag- 
nitude, the internal structure of the four populations in terms of stellar components is made 
more evident. In the i^-band we can clearly identify three peaks moving from right to left 
as -fCthr is increased from 20 to 30. As usual, the peak on the left is due to the tip of the 
RGB, the central peak is due to the red clump and the one on the right is made up of TOP 
stars, coming into play at higher thresholds. In the l^-band there is a first peak formed by 
red clump stars moving from high to low values of vg. At high thresholds the TOP stars are 
able to form a shoulder and then a second peak on the right of the red clump peak. At low 
values of Vthr the disk distribution is influenced by RSG, which form the tail at small vq. 
The contribution of RSG becomes subdominant already at intermediate thresholds, where 
it is taken over by red clump stars and finally by TOP stars. 



9. Rate and duration of the events 

An interesting quantity for the definition of the timescale of observational campaigns 
is the rate of the events, defined as the number of new detectable events occurring per unit 
time. It can be calculated by 

oo -Rz 27r 

Vp^ = j dMK np^ j dDi^s 2 J dr^ J dcj) fp^ ap^ G {vpjap^) , (21) 



where 

G{x) = y|e-^"/^ [(2 + x2)/o(xV4) + x'hix'/A)] (22) 

and is the modified Bessel function of order u. 

The combination crp^ G {vp^/crpj) is the result of the integration on the velocity distri- 
bution, again assumed to be a gaussian with dispersion a p . centered on the rotation velo city 



vp^. For further details on this point, see the appendix of [Alexander &: SternbergI (1l999l ). 
The function G{x) satisfies the following limits 

lim cTp, G ivpjap) = a/tt/^ ap (23) 

tip— >0 J ^ J J 

hm ap_ G (vpJc^p^ = vp (24) 
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Applying equation fl2T|) to each of our four populations, we obtain the curves in Figure 
[TH We see that the more distant populations move more slowly and thus generate new 
events at a lower rate. For this reason the rate of the bulge events practically coincides with 
the rate of P1+P2, while the rate of the disk drops to levels comparable to those of P3. 

Having calculated both the rate and the number of simultaneously present events, we 
immediately derive the average time spent by an event above threshold as 

< ATp, >= ^ ^ (25) 

where the last equality holds approximatively and involves the average radius of the lensing 
zone < Rz >, and the average transverse velocity < v± >. 

The average duration of the events for each source population are drawn in Figure [13 
It grows with the chosen threshold of the experiment. This is due to the fact that as Kthr 
is increased the lensing zone becomes larger and larger and the time required by the sources 
to cross it becomes proportionally larger (see the approximate expression in eq. [25]). 

The modulation in the growth at intermediate thresholds in the K-hand is again an 
effect of the presence of different stellar branches in each population. In fact, at low i^thr 
most of the events are due to bright sources in the red clump. Given their intrinsic luminosity, 
these sources generate secondary images that stay longer above threshold. At intermediate 
thresholds, there emerges a non-negligible contribution from fainter sources. These sources 
have a smaller lensing zone and their secondary images spend much less time above threshold 
with respect to the brighter sources. The distribution of the average duration becomes 
bimodal, with the brighter sources enjoying a long time above threshold and the fainter 
sources spending less time above threshold. The average time of such a bimodal distribution 
is therefore in the middle between the two maxima. As we increase iiTthr further, the fainter 
sources dominate with respect to the brighter and ATp^^ is determined by these sources only. 
Then the modulation observed at intermediate i^'thr in the curves in Figure [15] is just the 
transition from a regime dominated by more luminous sources to a regime dominated by 
dimmer sources. In the V-band this transition is less evident. 

Coming to a quantitative analysis of the average duration, we note that the events 
involving sources in the disk of M31 last up to several centuries. For the bulge the average 
time spent above threshold ranges from 50 to 200 years. For P1+P2 the situation is more 
dynamical, thanks to the higher velocities characterizing the stars very close to the central 
black hole. An observational campaign lasting several years could study the evolution of the 
secondary images already present and hope to see new events (if the threshold magnitude is 
sufficiently high). The average duration of the events involving sources in P3 is of the order 
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Fig. 14. — (a) Rate of the gravitational lensing events for the four source populations con- 
sidered in the text expected in the K-band as function of the threshold magnitude, (b) The 
same in the F-band. The rate of P1-I-P2 and the rate of the bulge practically coincide. 
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Fig. 15. — (a) Average time spent above threshold in the i^-band for the gravitational lensing 
events for the four source populations, (b) The same in the l^-band. 
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of one year, but the extremely low rate of new events for this population leads one to discard 
it as a source of gravitational lensing events. 



10. Identification of the lensing images 

On the basis of the analysis given in the previous sections, it is time to set up a realistic 
methodology to select the candidate lensing images and test their authenticity. 

As mentioned before, the most probable lensing events will involve bulge stars as sources, 
at typical distances ranging from a few pc to 1 kpc from the lens. In this configuration, the 
secondary images may form at an angular distance from the central black hole ranging from 
a few mas to 0.1". In order to get access to the bulk of the gravitational lensing events, it 
is then mandatory to employ interferometry. At this point, the question arises whether a 
secondary image of a lensed star can be distinguished from a background star unambiguously. 

In this respect, it is instructive to calculate the expected number of background stars 
in an angular area within radius 6 centered on the black hole. For a threshold magnitude 
Xthr, only the sources with observed magnitude K < Xthr will be observable. This fixes a 
threshold on the intrinsic magnitude of the sources as 

Mi^,max = -fi'thr " 5 logiQ - Ak- (26) 

The number of background sources is then 

Mk 



,(i^thr,^)= j dMKnp.(MK) j dz j dr r j d(l>fp.(z,r,(l>). (27) 



-oo --DoL 



Actually, since the spatial distributions for all source populations are centered on the 
black hole and have an extent much smaller than Dql, we can take the approximation 
z -C Dqi^ and push the lower extremum in the ^-integral to — oo. 

Let us imagine an observational campaign to discover secondary images due to grav- 
itational lensing lying within 6 — 0.01" of the central black hole. Fixing A'thr = 24, the 
number of background stars in such area is Bp^ = 3.5 x 10~^, i^pi+p2 = 8.3, i^Buige — 3-7, 

^Disk = 1.1 X 10-2. 

Our image will thus contain something like 12 stars within 0.01" from the central black 
hole. The first step would be to check for possible alignments between these stars and the 



- 39 - 



central black hole. If we find that the line joining two stars passes through the central 
black hole, we can select these stars as candidate lensing images of the same source. Since 
gravitational lensing is an achromatic phenomenon, the two images must have the same 
color indices (unless there is reason to believe that the photon paths of the two images cross 
regions with sensibly different amount of dust). Therefore, false candidates could be excluded 
taking images in different IR bands and checking whether the luminosity ratio between the 
two images remains the same in all bands. 

A further test is to take another image after several months or even the next year. If the 
alignment is just a chance product and one of the two stars has a high enough proper motion 
(which is likely to happen if the star is intrinsically close to the central black hole), then 
the alignment will be lost in the second image. Conversely, if the alignment is an authentic 
product of gravitational lensing, and if any evolution shows up, then the two images must 
move in such a way that the line joining them always passes through the black hole. However, 
as shown in in §[8], the proper motion of the secondary images is typically below 0.1 mas/y. 
Therefore, it would be improbable that the observed events show any evolution within less 
than ten years. Moreover, given the very low rate of new events, follow-up observations 
would make sense only if thresholds higher than i^thr = 30 are reached in the observations. 

Our estimate of 12 background stars in an area of radius 0.01" centered on the black 
hole is based on our spatial distributions, which do not take into account a possible presence 
of a further cuspy cluster of stars internal to P3 and thus very close to the supermassive 
black hole. If such a cluster were present, then the number of background stars would 
sensibly increase. However, these stars would move very fast and chance alignments in a 
first image would be easily discarded after the analysis of a second image. As stressed 
before, the most delicate background actually comes from stars that lie along the line of 
sight but intrinsically far from the black hole and consequently have slow proper motion. 
We believe that the number of these stars is correctly estimated through equation (1271) . For 
completeness, we also mention that the estimate of background stars for i^'thr = 30 yields 
1600 stars and for Vthr = 33 gives 1320 stars. 

As a final consideration, we note that the deformation of the point-like image of the 
source is of the order of 10~^ arcsec for the gravitational lensing events considered in this 
work. This prevents from using image deformation as a selection tool for lensing effects. 
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11. Conclusions 

In this paper we liave undertaken a deep investigation of the possible observation of 
gravitational lensing effects due to the black hole in the center of M31. We have considered 
stars belonging to four different populations (bulge, disk and the central clusters P1+P2 and 
P3) as candidate sources. Through detailed modelling of the spatial distributions and the 
luminosity functions of these four populations we have calculated several quantities that can 
be used to quantify the power of the central black hole as a gravitational lens using present 
and future observational facilities. We have carried out our analysis both in the i^-band and 
in the V^-band. 

The main outcome is the number of expected lensed sources at a given time whose 
secondary image is brighter than a specified threshold magnitude. The results are summa- 
rized in Figure [6] and Tables [3] and HI showing that indeed we expect 1.4 lensed sources at a 
threshold /^thr = 24, 16 events at i^'thr = 27 and 180 events at i^'thr = 30. In the \^-band we 
would have 1.3 events at Vftr = 27, 25 events at Vthr = 30 and 270 events at Vthr = 33. 

We have also presented the distribution of the lensed sources as a function of their 
absolute magnitude, showing the contribution of different stellar branches as sources of 
gravitational lensing events. 

The distribution of the lensed sources as a function of the distance shows that the disk 
stars are mostly lensed at about 1 kpc distance from the black hole, whereas the bulge stars 
can be lensed at distances ranging from a few pc to 1 kpc. The inner cluster stars are lensed 
at fractions of a pc. 

This difference is reflected in the angular separation of the images from the black hole, 
which ranges from 1 mas to 0.1" for bulge and disk stars, whereas it stays of the order 1 mas 
or below for the inner clusters. The resolution needed to get a significant number of events 
is thus of the order of a few mas, requiring the employment of long baseline interferometers 
or extremely large telescopes. This justifies our choice to present our analysis in the K- 
band, for which the present infrared interferometers such as Keck, LET and the future 
space telescope JWST are optimized. The VLTI, though representing the most advanced 
interferometer operating in IR bands, cannot observe M31 efficiently. The parallel analysis 
in the \^-band provides a deeper comprehension of all the effects coming into play. In this 
band, the expectations for the number of lensing events are comparable to those in the K- 
band, provided that V^-band interferometers are constructed in the future with sizes of the 
same order as in the i^-band. 

The gravitational lensing events discussed in this paper have a very slow evolution, 
with typical angular velocities of the images between 0.01 and 0.1 mas/y. We have outlined 
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a possible observational strategy to select gravitational lensing events within the expected 
background stars, by checking the alignment of pairs of images, the achromaticity of the flux 
ratio and setting up follow-up observations to detect any secular development. 

Our analysis can be repeated with different models of the source populations. Given 
the present uncertainties on the morphology of the components of M31, we believe that our 
estimates can be corrected by a factor of a few at most. 

The chances for concrete observations of gravitational lensing events by the supermassive 
black hole in M31 are rather low (though not null) even with the best facilities available today. 
The perspectives will be definitely increased with the realization of future projects such as 
the JWST, the extremely large ground telescopes or the realization of new long baseline 
optical interferometers in the northern hemisphere. With such facilities it should be possible 
to reach the higher thresholds indicated in this paper while keeping a very high angular 
resolution. This would open the way to an intensive research of gravitational lensing events 
and to their use in the investigation of the environment of the supermassive black hole. 

We are grateful to Gaetano Scarpctta and an anonymous referee for some fundamental 
comments which have taken to a consistent improvement of the manuscript. This work has 
made use of the lAG-STAR Synthetic GMD computation code. lAG-STAR is supported and 
maintained by the computer division of the Instituto de Astrofisica de Ganarias. The authors 
acknowledge support for this work by MIUR through PRIN 2006 Prot. 2006023491_003 and 
by research funds of the Salerno University. 
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